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ABSTRACT 

As part of the NASA EPOXI Mission of Opportunity, we observed seven known 
transiting extrasolar planet systems in order to construct time series photometry of 
extremely high phase coverage and precision. Here we present the results for four "hot- 
Jupiter systems" with near-solar stars— HAT-P-4, TrES-3, TrES-2 and WASP-3. We 
observe ten transits of HAT-P-4, estimating the planet radius R p = 1.332 ± 0.052 i?j up , 
the stellar radius it!* = 1.602 ± 0.061 Rq, the inclination i = 89.67 ± 0.30 degrees and 
the transit duration from first to fourth contact r = 255.6 ±1.9 minutes. For TrES-3, 
we observe seven transits, and find R p = 1.320 ± 0.057 i?j up , R* = 0.817 ± 0.022 R & , 
i = 81.99±0.30 degrees and r = 81.9±1.1 minutes. We also note a long term variability 
in the TrES-3 light curve, which may be due to star spots. We observe nine transits of 
TrES-2, and find R p = 1.169 ± 0.034 R Jnp , R± = 0.940 ± 0.026 R Q , i = 84.15 ± 0.16 
degrees and r = 107.3 ± 1.1 minutes. Finally we observe eight transits of WASP-3, 
finding R p = 1.385 ± 0.060 R Jup , R± = 1.354 ± 0.056 R Q , i = 84.22 ± 0.81 degrees and 
t = 167.3 db 1.3 minutes. We present refined orbital periods and times of transit for each 
target. We state 95% confidence upper limits on the secondary eclipse depths in our 
broadband visible bandpass centered on 650 nm. These limits are 0.073% for HAT-P-4, 
0.062% for TrES-3, 0.16% for TrES-2 and 0.11% for WASP-3. We combine the TrES-3 
secondary eclipse information with the existing published data and confirm that the 
atmosphere likely does not have a temperature inversion. 
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1. Introduction 



The EPOXI Mission of Opportunity is a re-purposing of the Deep Impact flyby spacecraft, and 
comprises the Extrasolar Planet Observation and Characterization (EPOCh) investigation and the 
Deep Impact extended Investigation (DIXI). The primary goal of EPOCh was to scrutinize a 
small set of known transiting extrasolar pla nets. From 2008 Janua ry to 2008 August, we used the 
high resolution imaging (HRI) instrument ([Hampton et al.l 120051 ) and a broadband visible filter 
to construct high precision, high phase coverage and high cadence light curves for seven targets. 
We observed each target nearly continuously for several weeks at a time. The main science goals 
of EPOCh were to refine the system parameters of the known planets, to search for additional 
planets both directly (via transits of the additional body) and indirectly (via induced changes in 
the transits of the known planet), and to constrain the reflected light from the known planet at 
secondary eclipse. It is also useful to provide updated periods and times of epoch for these systems 
in order to reduce uncertainties on predicted transit and eclipse times, and therefore maximize 
the return of follow-up observations. In previous EPOCh pa pers we have presented the search for 
additional pla nets in the GJ 436 system ([Ballard et al.ll2010l ) and the secondary eclipse constraints 
for HAT-P-7 ([Christiansen et al.ll2010l ). In this paper we present the updated system parameters, 
including constraints on the transit timing and changes in the transit parameters, and secondary 
eclipse constraints for a further four targets: HAT-P-4, TrES-3, TrES-2 and WASP-3, introduced 
below. The search for additional planets in these systems will be presented in a separate paper 
(Ballard et al. in prep). 



The exoplanet HAT-P-4b ([Kovacs et al.1120071 ) . orbits a slightly evolved metal-rich late F star. 
With a mass of 0.68 Mj up and a radius of 1.27 i?j up , it joined the ranks of inflated planets that 
have continued to challenge models of the physical structure of hot Jupiters. 



TrES-3 (jO'Donovan et al.l 120071 ) is notable for its very short orbital period of 1.30619 days. 
This proximity to the star makes TrES-3 a promising target for observations of reflected light at 
visible wavelengths; the planet-to-star flux ratio as measured in reflected light during the secondary 
eclipse is given by A g (R p /a) 2 , where A g is the ge ometric albedo, R„ i s the planetary radius and 



a is the semi- major axi s of the planetary orbit. IWinn et al.1 (|2008l ). Ide Mooii &: Snellen! ((2009) 
and iFressin et al.1 (j2010i ) have observed secondary eclipses of TrES-3 at visible and near-infrared 
wavelengths, and the emerging picture of the planetary atmosphere is one with efficient day-night 
re-circularization and no temperature inversion in the upper atmosphere. This is in contrast t o 
predictions of a tem perature inversion based on the high level of irradiation ([Fortnev et al.l 12008) . 



Sozzetti et al.l ([20090 studied the transit timing variations of TrES-3 and noted significant outliers 
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from a constant period. I Gibson et al.l (|2009l ) monitored further transit times of TrES-3 and ruled 
out sub-Earth mass planets in the exterior and interior 2:1 resonances for circular orbits. 



TrES-2 (jO'Donovan et alJl2006f) was the fir st transiting p l anet found in the field of view of 
the NASA Kepler mission (jBorucki et al.ll2009l ). iHolman et al.1 (|2007l ) noted that the high impact 
parameter (6 ~ 0.85) of TrES-2 made transit parameters such as inc lination and d urati on sensitive 
to changes due to orbital precession. iMislis Schmittl (|2009l ) and iMislis et al.1 (120101) claimed a 



signifi cantly shorter duration for TrES-2 transits two years after the measurements of lHolman et al 



(|2007l ). They proposed that this was caused by a change in orbit inclination due to precession, and 
that the durati on would continue t o dec rease. However. IScuderi et al.l (|2009l ) measured a duration 
consistent with lO 'Donovan et al.l (|2006l ) and IHolman et al.l (|2007l ) and did not see the predicted 
trend of decreasing transit dura tion. Secondary eclipses of TrES-2 have been observed in the near- 
infrared (jO'Donovan et al.ll2010l ). and the results favor a thermal inversion in the upper atmosphere, 
supporting the hypothesis that highly irradiated p lanetary atmospher es ha ve inversions. The t ransit 
timing variations of TrES-2 have been studied bv lRaetz et al.1 (|2009l ) and lRabus et al.1 (|2009l ). who 
find no statistically significant variations. 



WASP-3b (|Pollacco et al.l 12008 ). with a short period (1.84634 days) and a hot host star (F7- 
F8V, T e ff = 6400K), is one of the hottest transiting planets known, and another very good target 
for observing reflected light at secondary eclipse. 

The paper is organized as follows. The observations and generation of the light curves are 
described in Section [2j the transit analysis is presented in Section [3l the secondary eclipse analysis 
is presented in Section [J] and the results are discussed in Section 



2. Observations and analysis 

The EPOCh observations were made using the high resolution imager (HRI), which has a 30-cm 
aperture and a 1024x1024 pixel CCD. For our observations we used a clear visible filter, covering 
350-1000nm, in order to maximize the throughput of photons. The integration time for the science 
observations was 50 seconds, which for near-continuous observations results in roughly 1500 images 
per day. Since the on-board spacecraft memory is only 300Mb, we initially chose to read out only 
a 128 x 128 pixel sub-array of the full CCD, to ensure full phase coverage between data downlinks 
from the spacecraft. The CCD comprises four quadrants that are read out independently, and the 
sub-array is centered on the CCD where the four quadrants meet. The pixel scale is 0.4 arcsec 
pixel -1 , resulting in a sub-array field of view of 0.72 square arcminutes. The images are significant 
defocused, resulting in a stellar point-spread function (PSF) with a full-width half maximum of 4 
arcseconds. Typically this meant that the target star was the only star in the field of view, and 
we were unable to employ relative photometry techniques for removing correlated noise in the light 
curves. 

Table Q] summarizes the observing schedules for each of the four targets. HAT-P-4 and TrES-3, 
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along with GJ436 and XO-2, were observed during the initial observing block from 2008 January 
to 2008 May. The project was awarded an additional contingent observing block from 2008 June 
to 2008 August, during which time HAT-P-4 was re-observed, and TrES-2 and WASP-3 were also 
observed. During the contingent observations we began observing in a larger 256 x 256 pixel sub- 
array mode, to reduce losses from pointing drifts that occasionally resulted in the target star lying 
outside of the 128 x 128 pixel sub-array field of view. The number of images that could be obtained 
with the larger sub-array mode between data downlinks from the spacecraft was constrained by 
the data storage capacity on board the spacecraft. Therefore, in order to maximize the phase 
coverage we chose to restrict observations in the 256 x 256 pixel sub-array mode to the times of 
particular interest — during the transits and secondary eclipses. One event per data downlink could 
be observed in the larger sub-array mode without reducing the temporal coverage. Table [T] shows 
the total number of transits and eclipses observed for each target, with the number observed in the 
256 x 256 pixel sub-array mode given in parentheses. As discussed in Section [Tj TrES-2 was claimed 
to show changes in the transit inclination with time. Therefore, we used the larger sub-array mode 
to observe the transits of TrES-2 where possible. WASP-3 was a promising target for secondary 
eclipse observations, and therefore we observed the secondary eclipses of WASP-3 in the larger 
mode where possible. For HAT-P-4 we observed two of the three transits and two of the three 
eclipses obtained in the contingent observations in the 256 x 256 pixel sub-array mode. TrES-3 was 
observed in the initial observing block and no observations were obtained in the larger mode. 



2.1. Image Calibration and Time Series Extraction 



We receive calib rated FITS images from the extant Deep Impact data reduction pipeline 
(jKlaasen et al.l 120051 ) . These data have been bias- and dark-subtracted and flat-fielded, using 
calibration images obtained on the ground before launch. Due to the very high precision required 
in the light curves, we perform several additional calibration steps to account for changes in the 
CCD since launch. The spacecraft pointing drifts considerably with time, resulting in significant 
coverage of the CCD by the stellar PSF and plac ing paramount importance on the flat-fielding. 
The procedure is described in lBallard et al.l (|2010l ) and summarized here. 



Table 1. EPOCh observations 



Target 


VMag 


UT Dates observed (2008) No. 


of Transits 51 


No. of Eclipses a 


HAT-P-4 


11.22 


01/22-02/12, 06/29-07/07 


10 (2) 


9(2) 


TrES-3 


11.18 


03/06-03/18 


7(0) 


6(0) 


TrES-2 


11.41 


06/27-06/28, 07/19-07/29 


9 (7) 


8 (2) 


WASP-3 


10.64 


07/17-07/18, 07/30-08/07, 08/10-08/15 


8 (0) 


9 (8) 



a Including partial events. The number in brackets is the subset of events observed in 256 X 256 pixel 
sub-array mode. 
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For each target, we use a PSF constructed from the images to locate the position of the star to 
a hundredth of a pixel. At this stage we reject images with IOct outliers from the PSF fit, assuming 
the stellar PSF to be contaminated by an energetic particle hit. We subtract a time-dependent bias 
calculated for each quadrant from the corresponding overscan region. We reduce the pixels in the 
central columns and rows of the CCD (forming the internal boundaries between the quadrants) by 
roughly 15% and 1% respectively, to correct an artifact produced by the CCD readout electronics. 
For data obtained in the 256 x 256 pixel mode, we scale the images by a constant (typically differing 
from unity by one part in a thousand) to correct an observed flux offset between the two sub-array 
modes. 

In order to track time-dependent changes in the flat-field since launch, there is a small green 
LED stimulation lamp that can be switched on to illuminate the CCD. We obtained blocks of 200 
calibration frames using this lamp, which were taken every few days throughout the observations, 
alternating between blocks in the smaller and larger sub-array modes in the contingent observations. 
We correct each science frame by the flat-field generated from lamp images taken in the same sub- 
array mode. We assume any remaining flat-field errors to be color-dependent and therefore unable 
to be addressed by the monochromatic lamp. 

We perform aperture photometry, using a circular aperture of radius 10 pixels. The resulting 
light curves exhibit significant correlated noise on the order of 1%, which is associated with the 
drift in the spacecraft pointing. In order to correct for this, we use the data itself to generate a 
sensitivity map of the CCD. We assume the out-of-transit and out-of-eclipse data to be of uniform 
brightness, with two caveats. First, the star may have intrinsic variations in stellar brightness due 
to spots. Only one of the four targets displayed long-period variability (Figure [3]), and this was 
removed by fitting and removing a polynomial in time before producing the CCD sensitivity map. 
Second, transits of a dditional planets may be present, which will be suppressed with this treatment 



([Ballard et al.l |2010| ) . We randomly draw several thousand of the out-of-transit and out-of-eclipse 
points and find a robust average flux of the 30 spatially nearest neighbors. We use this set of 
averages to generate a two-dimensional surface spline to the flux distribution across the CCD. Each 
point in the light curve is then corrected by interpolating onto this surface. The entire procedure 
is iterated several times to converge on the positions and scaling factors that result in the lowest 
scatter in the out-of-transit and out-of-eclipse data in the final light curve. 

The robustness of the surface spline for each target depends on the coverage of the CCD by 
that target. If the coverage is small and the corresponding density of photometry apertures high, 
then there is a high probability that the same pixel will be returned to multiple times over the 
observations. Having flux measurements separated in time reduces the influence of stellar activity 
on our calibration of the sensitivity of each pixel. Figure [1] shows the complete CCD coverage for 
two targets. TrES-2 is well confined on the CCD and the density of photometry apertures leads 
to a more robust surface spline. The TrES-2 light curve prior to and post the application of the 
surface spline is shown in Figure 01 On the other hand, the photometry apertures for WASP-3 
sample a much larger area of the CCD, and in addition many of the observations obtained in the 
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256x256 pixel sub-array mode do not overlay the central 128x128 pixel sub-array. The resulting 
surface spline is therefore more sensitive to noise introduced by stellar activity or systematics that 
are not an artifact of the pointing jitter. The WASP-3 light curve prior to and post the application 
of the surface spline is shown in Figure The lower panel of Figure [5] shows how the noise in the 
final calibrated WASP-3 light curve bins down compared with the expectation for Gaussian noise, 
and the poor quality of the data is due to the low density of the CCD coverage for WASP-3. 




Fig. 1. — The CCD positions of the photometry apertures for two targets. Left: TrES-2 is confined 
to the center of the CCD and therefore the same pixels are sampled well for creating a robust surface 
spline. Right: WASP-3 is spread over a much larger fraction of the CCD, including large excursions 
out of the central 128x128 pixel sub-array when observations were obtained in the 256x256 pixel 
sub-array. This reduces the quality of the surface spline and results in a larger component of 
correlated noise in the WASP-3 light curve. 



2.2. Details for each target 

The final HAT-P-4 light curve is shown in Figure El HAT-P-4 was the first EPOCh target 
observed, initially for 22 days from 2008 January 22 to 2008 February 12, during the original 
EPOCh target schedule, and again for 8 days from 2008 June 29 to 2008 July 7 during the contingent 
observations. Of the 45,320 images obtained of HAT-P-4, 5434 were discarded due to the star being 
either out of the field of view or too close to the edge of the CCD to measure accurate photometry, 
1305 were discarded due to energetic particle hits, and 76 were discarded due to readout smear, 
for a final total of 38,505 acceptable images. All of the data obtained in the initial run are in the 
128 x 128 pixel sub-array mode. Of the contingent data, two of the three transits and two of the 
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three eclipses are in the larger 256 x 256 pixel sub-array mode, and the remaining data are in the 
smaller mode. The bottom panel of Figure [2] shows how the scatter in the final light curve scales 
down with increasing bin size — for Gaussian noise the expectation is the scatter will decrease as 
1/vJV, where N is the number of points in the bin. 




Fig. 2. — Upper panel: The full HAT-P-4 EPOXI light curve. The left panel shows the original 
run of seven consecutive transits. The right panel shows the three transits observed five months 
later during the EPOCh contingent observations. In each panel the lower curve is before the first 
application of the surface spline and the upper curve is the final calibrated light curve. The red 
data points were obtained in the larger 256 x 256 pixel sub-array mode. Lower panel: The scatter 
in the out-of-transit data with increasing bin size (diamonds) and comparing to the expectation 
for Gaussian noise (1/y/N, where N is the number of points in the bin, shown as the solid line 
normalized to the unbinned value of the scatter). The points do not follow the line, indicating 
residual correlated noise in the light curve. 

The final TrES-3 light curve is shown in Figure [3l TrES-3 was the second EPOCh target 
observed, for 12 days from 2008 March 6 to 2008 March 18. The gap in the light curve from 2-5 
days is due to a 'pre-look' for the subsequent EPOCh target, XO-2, which was performed in order 
to refine the pointing for that target. We obtained a total of 14,195 images of TrES-3, of which 
we discarded 1165 due to the star being out of or too close to the edge of the field of view, 1632 
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due to energetic particle hits and 127 due to readout smear, leaving 11,271 images. We obtained 
all of the TrES-3 data in the 128 x 128 pixel sub-array mode. After the initial application of the 
two-dimensional surface spline a long timescale, low amplitude variability was evident in the light 
curve. This can be seen in the lower light curve in Figure In order to remove this variability we 
bin the out-of-transit data by two hours and fit with a time-dependent fifth-order polynomial for 
the data occurring later than 4.0 days. We divide out this feature before iterating over the previous 
steps to produce the final light curve. The polynomial is plotted on the lower light curve, and the 
final light curve is shown as the upper curve in Figure [3j As with HAT-P-4, the bottom panel of 
Figure [3] shows the noise properties of the data. 
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Fig. 3. — Upper panel: The TrES-3 EPOXI light curve. The gap from 2.5-5 days is during the 
pre- look for a subsequent target. Seven transits of TrES-3 were observed in total. The lowest 
light curve is prior to the first application of the surface spline, the middle light curve is after the 
application of the spline but prior to the removal of the time-dependent polynomial, and the upper 
light curve is the final calibrated data set. Lower panel: See Figure [2] for explanation. In the case 
of TrES-3, where all data were obtained in the smaller sub-array mode and the total time span is 
relatively short, the scatter bins down close to the expectation for Gaussian noise. 



The final TrES-2 light curve is shown in Figure HI We observed TrES-2 during the contingent 
EPOCh observations, from 2008 July 7 to 2008 July 30, in addition to a pre-look for pointing on 
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2008 June 28 and 29. In total, we obtained 31,210 images of TrES-2, with 1979 discarded due to 
the star lying out of or too close to the edge of the field of view, 1427 discarded due to energetic 
particle hits and 80 discarded due to readout smear, for a total of 27,724 acceptable images. We 
observed nine transits in total, including seven in the 256 x 256 pixel sub-array mode. The lower 
panel of Figure S] shows that correlated noise remains in the final light curve. 

EPOCh TrES-2 Time Series 
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Fig. 4. — Upper panel: The TrES-2 EPOXI light curve. The data obtained from days 1-3 are the 
pre- look, for refinement of the spacecraft pointing. From days 3-11 the spacecraft was observing a 
different target before returning to TrES-2 with updated pointing parameters. The gap from days 
21-23 spans the pre-look for the subsequent target. Nine transits of TrES-2 were observed in total. 
The lower curve is prior to the first application of the surface spline and the upper curve is the final 
calibrated light curve. The red data points were obtained in the larger 256 x 256 pixel sub-array 
mode. Lower panel: See Figure [2] for explanation. 

The final EPOCh light curve for WASP-3 is shown in Figure We observed WASP-3 during 
the contingent observations, from 2008 July 29 to 2008 August 16, with a pre-look from 2008 July 
17 to 2008 July 19. We obtained 24,015 images of WASP-3, of which we discarded 4,182 due to 
the star being out of or too close to the edge of the field of view, 403 due to energetic particle 
hits, and 808 due to readout smear, leaving 18,622 acceptable images. For WASP-3, none of the 
eight transits were observed in 256 x 256 pixel sub-array mode, however eight of the nine secondary 
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eclipses were observed in this mode. The two-dimensional surface spline relies on multiple visits to 
the same part of the CCD to characterize robustly the interpixel variations. This is particularly 
true for the data that occur during the transits and eclipses, since they cannot be assumed to be 
of uniform flux and are therefore excluded from the creation of the surface. In order to effectively 
flat-field the data that are taken during transit and eclipse, the observations taken during these 
times must be gathered at the same spatial positions as data obtained at other times. In the case of 
WASP-3, four of the eight secondary eclipses occurred at locations that were poorly sampled. No 
out-of-transit or out-of-eclipse observations fell on these pixels, and therefore we cannot estimate 
the true sensitivity of these pixels in order to produce an effective flat-field. These eclipses occur 
at 1.0, 17.6, 19.3 and 26.6 days, and can be seen in the light curve as increases in flux. These 
four eclipses are discarded for the final analysis. Besides these events, a significant fraction of the 
WASP-3 data are distributed in poorly-sampled areas of the CCD, degrading the robustness of the 
two-dimensional surface spline. The bottom panel of Figure demonstrates the adverse effect this 
has on the noise properties of the final light curve, as the data do not bin down as expected for 
Gaussian noise. 



3. Transit analysis 



For the transit analysis, we make several additional calibration steps. The two-dimensional 
surface spline uses only a fraction of the data to generate the surface, in order to preserve as much 
of the information in the light curve as possible, and to minimize the suppression of transits of 
putative additional planets. However for the transit analysis, we use all of the available data to 
calibrate each event. For each transit, we define a window approximately three times the duration 
of the transit, centered on the predicted transit time. We take each point in this window and divide 
the flux by a robust average of the 30 spatially nearest points that do not fall in any of the transit 
windows. This is essentially a point-by-point application of the full two-dimensional surface spline. 
We then fit a slope, linear with time, to the out-of-transit data across each transit and divide it 
out, to remove any residual long timescale trends. 

For TrES -3, TrES-2 an d WASP-3, we generate non-linear limb-darkening coefficients of the 
form given by IClaretl ( 200d ) . I^/Ii = 1 — Ylti=i c «(^ ~~ A 4 ™^ 2 ); where I\ is the specific intensity at 
the center of the disk and ix = cos(j), with 7 the angle between the emergent intensit y and the line 
of sight. We use photon- weighted stellar atmosphere models of iKuruczl (|1994l . 120051 ) that bracket 
the published values of stellar T e g and log g, and convolve these with the total EPOXI response 
function, including filter, optics and CCD response. We fit for the four coefficients of the non-linear 
form of the limb-darkening using 17 positions across the stellar limb, at 2 nm intervals along the 
350-1000 nm bandpass. We calculate the final set of coefficients as the weighted average when 
integrated over the bandpass, and bi- linearly interpolate across T e g and log g for each target. The 
final set of coefficient for each targets is given in Table [3] for TrES-3, Table U] for TrES-2 and Table 
EJfor WASP-3. 
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Fig. 5. — Upper panel: The WASP-3 EPOXI light curve. The first two days of data are the pre-look 
to refine the pointing. The gap between 22 and 24 days is due to the pre-look for the subsequent 
target. The significant positive deviations seen at 1.0, 17.6, 19.3 and 26.6 days are instrumental in 
nature; see the text for details. The lower curve is prior to the first application of the spatial spline 
and the upper curve is the final calibrated light curve. The red data points were obtained in the 
larger 256 x 256 pixel sub-array mode. Lower panel: See Figure [2] for explanation. 



The quality of the EPOCh light curves is nearly sufficient to fit for the limb-darkening coeffi- 
cients rather than assuming theoretical values. Ultimately, the degeneracies between the geometric 
parameters of the transiting system and the limb-darkening coefficients prevent us from placing 
meaningful constraints on the coefficients. In the case of HAT-P-4 however, the system is very close 
to edge-on (i = 89.9l££ de grees), which reduces the parameter space considerably. Therefore, for 
HAT-P-4 we instead use a quadratic equation for the limb-darkening, I^/Ii = 1 — ail — fj,) — 6(1— /x) 2 , 
and allow two linear combinations of the coefficients, c\ = 2a + b and C2 = a — 2b to be free param- 
eters in the transit analysis, which produced a better fit to the data as defined below. 



When fitting the transits, we use the analytic equations of lMandel fc Agoll ([2002 ) to generate 
a model transit, and use \ 2 as a goodness-of-fit estimator. We use the Levenberg-Marquardt 
algorithm to fit three dimensionless geometric parameters of the system: R p /R+, R*/a and cos i, 
where R p is the planetary radius, R* is the stellar radius, a is the semi-major axis of the planetary 
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orbit and % is the inclination of the orbit. We fix the period to the published value, but allow the time 
of center of transit to vary independently for each of the transits. We then use the published mass 
values for e ach of the sys t ems t o convert the tr ansit paramete r s to p hysical prop erties, drawing 
values from iKovacs et all (|2007h fo r HAT- P-4, ISozzetti et all (|2009h for TrES-3, ISozzetti et al 



(|2007|) for TrES-2 and IPollacco et al.1 ((2003) for WASP-3. The final results of these fits are given 
in Table [2] for HAT-P-4, Table E] for TrES-3, Table 1 for TrES-2 and Table [5] for WASP-3. We 
also give the transit duration from first to fourth contact for each best-fit model. For WASP-3, we 
discard the final transit (which was significantly offset in flux due to correlated noise), and also a 
partial transit (which included only the ingress), for a total of six transits. The phase- folded and 
binned transits for each target are shown in Figure [6] for HAT-P-4, Figure [7] for TrES-3, Figure [8] 
for TrES-2 and Figured for WASP-3. 

The errors on the pa rameters are calculated using the residual permutation "rosary bead" 
method (|Winn et al.ll2008l ). For each target, we find the residuals to the best-fit model. We shift 
these residuals forward collectively to the next time stamp and add the best fit models back to the 
new residuals, generating a new realization of the light curve which retains the correlated noise 
signals in the original light curve. We repeat this process 8000 times (covering approximately six 
days) and each time we fit for and record the geometric parameters, times of center of transit, and 
limb-darkening coefficients where appropriate. For each parameter we construct a histogram of the 
8000 measurements, to which we fit a Gaussian. We then define the error on that parameter by 
the half-width half-maximum value of the best fit Gaussian. We find that increasing the number 
of iterations beyond 4000 does not significantly change the calculated errors. 

To find the errors in the transit times, we perform a second rosary bead analysis, holding 
the geometric and limb-darkening values fixed and allowing only the times of center of transit to 
vary. We find that 4000 iterations are sufficient to sample the range of correlated noise signals, and 
calculate the errors in the same fashion as the geometric parameters. For each target we calculate 
a new orbital period and epoch by performing a weighted linear fit to the EPOCh transit times 
and any published transit times. 



4. Secondary eclipse constraints 

Our constraints of the secondary eclipse depths are limited by the correlated noise in the 
EPOXI data. Ideally, for each target we would combine our multiple observations of the secondary 
eclipses to increase the signal to noise. However, the fluctuations due to the correlated noise 
preclude this. For example, Figure [11] shows six of the TrES-3 secondary eclipses, where in some 
cases correlated noise results in an increase in flux at the time of secondary eclipse, instead of 
the expected decrement. If we assume that the secondary eclipse in the EPOCh bandpass, with a 
central wavelength of 650 nm, is due exclusively to the reflected light of the planet, then the eclipse 
depths we would anticipate, for a geometric albedo of 1, would range from 0.02% for HAT-P-4 to 
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Table 2. HAT-P-4 system parameters 



Parameter Value 



Adopted values' 1 

Ah (M ) 1.26 ±0.14 

M p (M Jup ) 0.68 ±0.04 



Transit fit values 
Rp/ R* 
a/R* 
i (deg) 

Derived values 
P (days) 
T c (BJD) 
R* {Re) 
Rp {Rjup) 

t (mins) 

Limb-darkening coefficients 

a 

b 

Transit times (BJD) 



0.0855 ± 0.0078 
0.1672 ±0.0078 
89.67 ±0.30 



3.0565114 ±0.0000028 
2, 454, 502.56227 ± 0.00021 
1.602 ±0.061 
1.332 ±0.052 
255.6 ± 1.9 



0.314 
0.366 

2, 454, 490.33445 ± 0.00072 
2, 454, 493.39232 ± 0.00061 
2, 454, 496.44984 ± 0.00056 
2, 454, 499.50426 ± 0.00070 
2, 454, 502.56156 ± 0.00056 
2, 454, 505.62006 ± 0.00082 
2, 454, 508.67569 ± 0.00056 
2, 454, 649.27624 ± 0.00064 
2, 454, 652.33053 ± 0.00065 
2, 454, 655.38842 ± 0.00065 



a Masses are from lKovacs et al.l j2007l l. 
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Table 3. TrES-3 system parameters 



Parameter 



Value 



Adopted values a 
Ah (M ) 
M p (M Jup ) 

Transit fit values 
Rp/ Ftk 
a/K* 
i (deg) 



n Q98 +0.028 
0.928_ Q 04g 

1 910+ ' 075 



0.1661 ±0.0343 
0.1664 ±0.0204 
81.99 ±0.30 



Derived values 
P (days) 
T c (BJD) 
R* (Re) 
Rp (Rjup) 

i (deg) 
t (mins) 



1.30618608 ±0.00000038 
2, 454, 538.58069 ± 0.00021 
0.817 ±0.022 
1.320 ±0.057 
81.99 ±0.30 
81.9 ± 1.1 



Limb-darkening coefficients 

Cl 
C2 
C3 
C4 



0.5169 
-0.6008 
1.4646 
-0.5743 



Transit times (BJD) 



454, 532.04939 
454,533.35515 
454, 537.27463 
454,538.58126 
454, 539.88703 
454,541.19261 
454, 542.49930 



± 0.00033 
± 0.00035 
± 0.00038 
± 0.00035 
± 0.00040 
± 0.00035 
± 0.00041 



a Masses are from lSozzetti et al. 
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Table 4. TrES-2 system parameters 



Parameter Value 



Adopted values a 

hU (M Q ) 0.98 ±0.062 

M p (Afj U p) 1.198 ±0.053 



Transit fit values 
Rpl R* 
a/R* 
i (deg) 

Derived values 
P (days) 
T c (BJD) 
R* (Rq) 
Rp (Rjup) 

t (mins) 

Limb-darkening coefBcients 

ci 

C2 

C3 

C4 

Transit times (BJD) 



0.1278 ±0.0094 
0.1230 ±0.0179 
84.15 ±0.16 



2.47061344 ± 0.0000075 
2, 454, 4664.23039 ± 0.00018 
0.940 ± 0.026 
1.169 ±0.034 
107.3 ± 1.1 



0.3899 
-0.1391 

0.9662 
-0.0329 

2, 454, 646.93735 ± 0.00032 
2, 454, 656.81879 ± 0.00034 
2, 454, 659.28871 ± 0.00042 
2, 454, 661.76005 ± 0.00044 
2, 454, 664.23072 ± 0.00050 
2,454, 669.17156 ±0.00028 
2, 454, 671.64117 ± 0.00028 
2,454, 674.11318 ±0.00033 
2,454, 676.58257 ±0.00051 



a Masses are from lSozzetti et al.l ll2007h . 
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Table 5. WASP-3 system parameters 



Parameter 



Value 



Adopted values a 

(Mq) 
M p (M Jup ) 



1.24 
1.76 



+0.06 
0.11 
+0.08 
0.14 



Transit fit values 

Rp/Ri, 

a/R* 

i (deg) 



0.1051 ±0.0124 
0.1989 ±0.0287 
84.15 ±0.16 



Derived values 
P (days) 
T c (BJD) 
R* (Re) 
Rp (Rjup) 

i (deg) 
t (mins) 



1.8468373 ±0.0000014 
2, 454, 686.82069 ± 0.00039 
1.354 ±0.056 
1.385 ±0.060 
84.22 ±0.81 
167.3 ± 1.3 



Limb-darkening 

ci 
C2 
C3 
C4 



coefficients 



0.2185 
0.6183 
-0.1040 
-0.0426 



Transit times (BJD) 



454, 
454, 
454, 
454, 
454, 
454, 
454, 
454, 



679.43264 
681.27911 
683.12740 
684.97486 
686.82053 
690.51381 
692.36117 
694.20711 



± 0.00050 
± 0.00040 
± 0.00035 
± 0.00027 
± 0.00059 
± 0.00055 
± 0.00043 
± 0.00042 



a Masses are from lPollacco et al. I tosh . 
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Fig. 6. — The seven HAT-P-4 transits from the original observing schedule, phase-folded and binned 
in five minute intervals, light curve. The solid line is the best fit transit model. Lower panel: The 
residuals when the best-fit model is subtracted from the data. 

0.08% for TrES-30 

Since the fluctuations from correlated noise in the measured eclipse depths are sometimes 
larger than the signal we expect to measure, we choose not to combine the multiple observations 
and instead analyze each eclipse independently. Our intent is to use the scatter of individual eclipse 
measurements to constrain the amplitude of the correlated noise. As for the transits, for each eclipse 
in the data we apply a point-by-point correction to the data in and adjacent to the eclipse. For the 
targets presented here we assume that e = and therefore that the secondary eclipse occurs at a 
phase of 0.5. For TrES-3 and TrES-2 this assumption is strongly supported by previous secondary 
eclipse measurements with the Spitzer IRAC instrument which demonstrated no evidence of non- 
zero eccentricity (Fressin et al. 2010 and O'Donovan et al. 2010 respectively). For HAT-P-4 and 
WASP-3, the extant radial velocity data are consistent with circular orbits (Kovacs et al. 2007 



In fact, the CCD is found to be quite efficient at the redder wavelengths, and it is therefore feasible that for the 
hottest planets there may be a contribution from the thermal emission of the planet, resulting in deeper secondary 
eclipses. 



- 18 - 




-0.004 F ^ 

-0.06 -0.04 -0.02 0.00 0.02 0.04 0.06 

Orbital Phase 

Fig. 7. — Upper panel: The seven TrES-3 transits, phase-folded and binned in two minute intervals. 
The solid line is the best-fit transit model. Lower panel: The residuals when the best-fit model is 
subtracted from the data. 

and Pollacco et al. 2008 respectively). In addition to the point-by-point correction, we fit a linear 
time-dependent slope to the adjacent out-of-eclipse data to remove any remaining long timescale 
trends. Finally, we separate the data into 10 minute bins and remove 3a flux outliers from each 
bin. For TrES-2, we discard the first observed eclipse, which was obtained during the pre-look for 
this target, since the pre-look data are not well calibrated by the surface spline generated for the 
remaining data. We assume this is due to changes in the CCD in the time that occurred between 
the pre-look and the full set of observations. We also discard eclipses where less than half the event 
is observed, one for TrES-2, one for WASP-3 and one for HAT-P-4. As discussed in Section I2.2[ 
we finally discard four of the nine WASP-3 secondary eclipses that fall on regions of the CCD we 
cannot calibrate. 

We fit the eclipses using a transit model with the best-fit parameters from the transit analysis 
and no limb-darkening. We then scale the depth of this model to fit the data, finding the depth 
that minimizes the x 2 value. For each target we then find the mean (x) and standard deviation 
(a x ) of the individual best-fit depths, and define the 95% confidence upper limit on the eclipse 
depth as x + 2a x . The secondary eclipses of HAT-P-4, TrES-3, TrES-2 and WASP-3 are shown in 
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Fig. 8. — Upper panel: The nine TrES-2 transits, phase-folded and binned in two minute intervals. 
The solid line is the best-fit transit model. Lower panel: The residuals when the best-fit model is 
subtracted from the data. 



Figures 9-12. The upper limits are given in Table We note that we achieve a useful constraint 
only in the case of TrES-3. 



Discussion 



5.1. HAT-P-4 



For HAT-P-4, our estimates of the system par ameters are consisten t wit h (and in the case o f 
inclination, more precise than) those published by iKovacs et al.l (|2007l ) and I Torres et al.l ((2008). 



We calculate R p = 1.332 ± 0.052 R Jup , R± = 1.602 ± 0.061 R Q , i = 89.67 ± 0.30 degrees and 
r = 255.6 ± 1.9 minutes, where r is the transit duration from first to fourth contact. We use the 
discovery epoch and the ten EPOCh transit times presented in this paper to produce a new refined 
ephemeris of T C (BJD) = 2454245.81531 ± 0.00021 + 3.0565114 ± 0.0000028£. Figure [H shows the 
residuals to the new ephemeris. We see no evidence for transit timing variations in the residuals 
which have a scatter of roughly 2 minutes. 
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Fig. 9. — Upper panel: The eight WASP-3 transits, phase-folded and binned in two minute intervals. 
The solid line is the best-fit transit model. Lower panel: The residuals when the best-fit model 
is subtracted from the data. The significant in-transit deviation from the model is discussed in 
Section [5j 



We use eight of the nine observed secondary eclipses to constrain the depth of the eclipse in 
the EPOCh bandpass, discarding the ninth due to poor coverage of the event. The eclipses are 
shown in Figure fTUl We set a 95% confidence upper limit on the eclipse depth of 0.073%, which, if 
it were produced entirely by reflected light, would correspond to a planetary geometric albedo of 
A g = 3.5, a physically impossible value. In the future, full phase curves of HAT-P-4 are scheduled 
to be observed in the near-infrared 3.6 and 4.5 micron IRAC bands, as part of the Warm Spitzer 
census of exoplanet atmospheres, at which point we may begin to study the atmosphere in more 
detail. 



5.2. TrES-3 



F o r TrES-3, we find sy stem parameters consiste nt with those published by I O' Donovan et al 
(|2007h . lSozzetti et al.l feoOfll ) and lGibson et al.l feoOflh . with R p = 1.320±0.057 i?j up , R* = 0.817 ± 
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Fig. 10. — Eight EPOCh secondary eclipse observations of HAT-P-4, offset in relative flux for 
clarity and binned in five minute intervals. The error on each point is a/y/N, where a is the scatter 
in the bin and N the number of points. The solid lines are the best fit eclipse model in each case. 
The bottom three eclipses were obtained in the contingent block of observations. 



0.022 Rq, i = 81.99 ± 0.30 degrees and r = 81.9 ± 1.1 minutes. In the upper panel of Figure [TBI 
we plot the published values of inclination with time and note a weak trend towards decreasing 
inclination, however it is present at only the 1.5a level, and hence not significant (and largely 
dependent on the most recent value from Sozzetti et al. 2009). A more model- independent way 
of constraining changes in the transit parameters with time is by measuring the transit duration. 
Where available, we use the quoted transit duration and error, and otherwise we calculate the 



Table 6. EPOCh secondary eclipse measurements 



Target Eclipse Depth Upper limit Implied A g 



HAT-P-4 -0.0069 ± 0.0397% 0.073% 3.5 

TrES-3 -0.020 ±0.041% 0.062% 0.81 

TrES-2 0.023 ±0.071% 0.16% 6.6 

WASP-3 0.023 ±0.044% 0.11% 2.5 
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Fig. 11. — Six EPOCh secondary eclipse observations of TrES-3, offset in relative flux for clarity. 
The error on each point is a/y/N, where a is the scatter in the bin and N the number of points. 
The solid lines are the best fit eclipse model in each case. 



transit duration from the published para meters, using equatio n (4) from lCharbonneau et al.l ([2006). 
Following the analytic approximation of ICarter et al.l ( 20081) . we set the error on these calculated 
transit durations to twice the error in the measured transit times for each source. Although this 
error was originally derived for the transit duration from mid- ingress to mid-egress, as compared 
to the transit duration from first to fourth contact, we find that for the EPOCh data the errors 
calculated using this approximation and the errors measured from the data themselves are nearly 
identical (1.1 minutes and 1.0 minutes respectively). We plot the derived values in the lower panel 
of Figure [151 an d we see no evidence of a change in the transit duration with time. 



In Section [221 w e noted that in the process of calibrating the light curve, a long term variability 
was evid ent. This variability is consistent with stellar variability due to spots. Using a usini of <2 
km s _1 ( O'Donovan et al.l 120071 ). the rotational period of TrES-3 must be >21 days, considerably 
longer than our observation span of 12 days. We can therefore not place any additional constraints 
on the rotational period of TrES-3, however we note that if the variability is due to spots on the 
stellar surface rotating in and out of view then additional monitoring of TrES-3 may reveal the 
rotational period. 



- 23 - 



1.00 - 



0.99 



0.98 - 



0.97 



0.96 



— i — i — i — | — i — i — i — | — i — i — i — | — i — i — i — | — i — i — i — | — i — r 
TrES-2, binned x 5 minutes 

^V+n V++ y+/V^ ' 



♦ ♦ ♦ ♦ 

S^WrrtV^ 

m "'.W^ ^ 



0.44 0.46 



0.48 0.50 0.52 
Phase 



0.54 



0.56 



Fig. 12. — Six EPOCh secondary eclipse observations of TrES-2, offset in relative flux for clarity. 
The error on each point is a/y/N, where a is the scatter in the bin and N the number of points. 
The solid lines are the best fit eclipse model in each case. 



For TrES-3, we calculate a new ephemeris of T C (BJD) = 2454538. 58069±0. 00021+1. 30618606± 
0. 00000038.E using the published transit times and the seven EPOCh transits presented in this 
paper. Figure [16] shows the residuals to the new ephemeris. We see no evidence of the period 
changing with time or transit timing variations larger than 1 minute. 

Using the six EPOCh secondary eclipse observations of TrES-3, shown in Figure [11] we set a 
95% confidence upper limit on the eclipse depth of 0.062%. Thi s indic ates the planetary geometric 
albedo must be A g < 0.81 in the EPOCh bandpass. I Winn et al.l (|2008l ) observed secondary eclipses 
of TrES-3 in the i, z and R bands, and were able to put 99% confidence upper limits on the 
eclipse depths of 0.024%, 0.050% and 0. 086% respectively. The EP OCh upper limit at 0.65 fj,m is 
consistent with the R band upper limit. Ide Mooii fc S nellen (2009 ) ob served the secondary eclipse 



Fressin et al 



(|2010i ) observed secondary eclipses 



in the K band and found a depth of 0.241±0.043%. 
of TrES-3 with the Spitzer IRAC instrument, measuring depths of 0.356±0.036%, 0.372±0.054%, 
0.449±0.097% and 0.475±0.046% in the 3.6, 4.5, 5.8 and 8.0 micron bands respectively. The 
secondary eclipse measurements are shown in Figure [T7] as a function of wavelength. 
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Fig. 13. — Four EPOCh secondary eclipse observations of WASP-3, offset in relative flux for clarity. 
The error on each point is a/y/N, where a is the scatter in the bin and N the number of points. 
The solid lines are the best fit eclipse model in each case. 



Given the high l evels of stellar irradiation, the atmosphere of TrES-3 was anticipated to host 



a thermal inversion ( Fortnev et al 



2008 



de Mooii h Snellen! 120091 ) . Using all data sets, however, 



Fressin et al.l ( 20101 ) found the observations to be best fit with a dayside atmosphere model without 
a thermal inversion. 



Our model spectra are computed using the exoplanet atmosphere model developed in lMadhusudhan h Seae 



(|2009l ) . The model consists of a line- by- line radiative transfer model, with constraints of hydrostatic 
equilibrium and global energy balance, and coupled to a parametric pressure-temperature (P-T) 
structure and parametric molecular abundances (parametrized as deviations from thermochemical 
equilibrium and solar abundances). Our modeling approach allows one to compute large ensembles 
of models, and efficiently explore the parameter space of molecular compositions and temperature 
structure. 

We confirm previous findings that existing detections of day-side observations can be explained 
to within the la uncertainties by models without thermal inversions. The black curve in Figure H71 
shows one such model spectrum, which has a chemical composition at thermochemical equilibrium 
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Fig. 14. — The transit times of HAT-P-4. The o pen diamonds are t he EPOCh transit times from 
this paper; the asterisk is the discovery epoch (jKovacs et al.l 120071 ). Lower panel: An expanded 
view of the EPOCh transit times, with la errors of 48-71 seconds. 



and solar abundances for the elements. The mo del is also consisten t with the EPOCh upper-limit 
at 0.65 microns, and with the upper-limits from I Winn et al.l (|2008l ). The dark green dashed curve 
shows a 1600K blackbody spectrum of the planet, indicating that the data cannot be explained 
by a pure blackbody. The model reported here has a day-night energy redistribution fraction of 
0.4, indicating very efficient redistribution. Therefore, based on previous studies and our current 
finding, existing data do not require the presence of a thermal inversion in TrES-3. However, a 
detailed exploration of the model p arameter space would be neede d to rule out thermal inversions 
with a given statistical significance (jMadhusudhan fe; Seagerll2010l ). 



5.3. TrES-2 



For TrES -2, we derive syst e m par a meters that are consis tent a t the 1.5a leve l with estimates 
published by lo 'Donovan et al.l ((200Q), ISozzetti et all feooj ). and iHolman et all \200i ). finding 
R p = 1.169 ± 0.034 R Jup , = 0.940 ± 0.026 R Q , i = 84.15 ± 0.16 degrees and r = 107.3 ± 1.1 
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Fig. 15. — Upper panel: The estimates of the inclination for TrES-3 as a function of time. Lower 
panel: The estimates of the TrES-3 transit durations. 



minutes. 



As discussed in Section [JJ there is currently a debate as to whether the inclination of the 
planetary orbit and duration of the TrES-2 transit are decreasing with time due to orbital precession. 
In the upper panel of Figure [TBI w e plot the estimates f o r the inclinat i on as a fun ction of time. For the 
inclination, the error bars of lMislis fc SchmittJ (fcoOflh . iMislis et all (hoioh and lScuderi et all fcoOflh 
were calculated by fixing the stellar and planetary radii and allowing only the inclination and time 
of center of transit to vary. The remainder of the inclination error bars were calculated allowing 
all of the geometric parameters to vary simultaneously, which explains why they are considerably 
larger than the later results. Since the errors skew any weighted linear fit towards an unrealistically 
large increase in the inclination with time, we instead plot an unweighted linear fit to guide the 
eye. We note that TrES-2 is in the Kepler field and that any change in inclination with time will 
soon be measured with exquisite precision. 

The inclination measured from a particular transit light curve will necessarily depend on the 
geometric parameters and to some extent the choice of limb-darkening treatment. However, the 
transit duration is directly measurable from the light curve and should not depend on the limb 
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Fig. 16 — T he times of transit o f TrES-3. O'D07: lO'Donovan et al.1 ((20071 ): S09: ISozzetti et al 



(|2009l ): G09: iGibson et al.l (|2009l ): EPOXI: this paper. Lower panel: An expanded view of the 
EPOCh transit times, with la errors of 28-35 seconds. 



darkening. The lower panel of Figure [T5] shows the published transit durations as a function of 
time. Where they were not given, we calculated the durations and errors as described for TrES-3. 
In this case, we perform a weighted linear fit and do see a formally significant decrease in the 
transit duration with time. How ever, this conc l usion is heavily dependent on one point, in this 
case the duration calculated from lHolman et al.1 (|2007l ). If this point is excluded from the fit, then 
dr/dt = — 0.0015±0.0015, consistent with no change in the transit duration with time and therefore 
we do not claim to have detected a change in the transit duration with time. Again, we expect 
Kepler to provide a clear answer to this question. 

Using the published transit times of TrES-2 and the nine transits observed by EPOCh pre- 
sented in this paper, we find a new weighted ephemeris of T C (BJD) = 24544664.23039 ± 0.00018 + 
2.47061344 ± 0.00000075^. The residuals to this ephemeris are shown in Figure [T9j In the EPOCh 
residuals, we see no variations in the transit times above the level of 2 minutes; excluding the 
amateur data from the Exoplanet Transit Database due to the large error bars, the scatter in the 
full set of residuals is less than 5 minutes. We see no evidence for long term drifts in the period. 
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Fig. 17. — The optical and near-infrared secondary eclipse measurements of TrES-3. The EPOCh 
upper li mit of 0.062% is sho wn in blue at 0.65 microns. The remainin g upper limits in the optica l 
are from I Winn et al.l (|2008l ): the measurement at 2.2 micro ns is from Ide Mooii Snellen! ((2009); 



and the four measurements from 3.6 to 8.0 microns are from Fressin et al 



1(201(1 ) . The solid black 



line is a representative model from the set of models that fit the data to within la, and the dashed 
lined shows a black-body spectrum for a temperature of 1600K. The green circles represent the 
model integrated to the Spitzer bandpasses. The inset is the temperature-pressure profile for the 
model shown. 



We used six of the eight EPOCh secondary eclipses of TrES-2 to place a 95% confidence upper 
limit on the eclipse depth of 0.16%. This corresponds to a planetary geometric albedo of A g = 6.6. 
As for HAT-P-4, this is not a physically plausible value. 



5.4. WASP-3 



For WASP-3, we measure system para meters that ar e cons isten t with, and an impr ovement 
upon, previously published parameters from lPollacco et al.l (|2008l ) and lGibson et al.l (|2008l ). finding 
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Scuderi et al. 



20091 ); EPOXI: this paper. Lower panel: The TrES-2 transit durations. In this case the dotted line 



is a weighted linear fit. 



R p = 1.385±0.060 Rjup, R* = 1.354±0.056 R Q , i = 84.22±0.81 degrees andr = 167.3±1.3 minutes. 
We generate a new refined ephemeris from the published transit times and the eight EPOCh transits 
in this paper, finding T C (BJD) = 2454686.82069 ±0.00039 + 1.8468373 ±0.0000014^. The residuals 
to this ephemeris are shown in Figure [20j 

The phase- folded light curve of WASP-3 (Figure [9]) shows correlated residuals in the latter half 
of transit. Since the noise in the transit exceeds the noise out of transit, one conclusion could be 
spot activity on the surface of the star being eclipsed during transit. However, if we examine the 
transits individually we observe that the correlated noise in the full light curve is not typically larger 
in transit than out of transit. The six transits used in the analysis are shown in Figure I2T1 In the 
transits numbered 2, 3, 4, and 6 large deviations can be seen in the second half of the transit, which 
leads to residuals in the phased light curve. If there were star spots producing correlated residuals 
in the transits, we would not necessarily expec t them to occu r at th e same phase for each transit. 



The v sin i for WASP-3 has been measured bv lSimpson et al.l (|2009l ) to be 15.7l}' 3 km s 1 , which 



- 30 - 



o 
l 

o 



20 P 



10 r 



-20 § 




I ~ I 

□ O'D06 xH07 ♦ R09 * R09(ETD) o EPOXI 



4000 



4200 4400 

BJD-2450000.0 



4600 



4800 



o 
I 

o 




4650 



4660 4665 
BJD-2450000.0 



4670 



4675 



4680 
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Holman et al.l ([20071); R09: iRaetz et al.l (|2009h : ROQfETD^ lRaetz et al.1 (|2009h (from the Exoplanet 
Transit Database, |http://var.astro.cz/ETD[ ); EPOXI: this paper. Lower panel: An expanded view 
of the EPOCh transit times, with la errors of 24-44 seconds. 



corresponds to a rotational period for the star of 4.2 days. Transits of WASP-3 are spaced by 1.85 
days, so it is improbable for spot activity to appear at the same phase in successive transits. Given 
these constraints, we conclude that the alignment of signals with phase in the EPOCh transits of 
WASP-3 are coincidental and are due to instrumental artifacts. 

We use four of the nine EPOCh secondary eclipse observations of WASP-3 to set a 95% 
confidence upper limit on the eclipse depth of 0.11%. This corresponds to a planetary geometric 
albedo of A g = 2.5 in the EPOCh bandpass, which is not a useful constraint for the planetary 
atmosphere. 



6. Conclusion 



We have presented time series photometry from the NASA EPOXI Mission of Opportunity 
for four known transiting planet systems: HAT-P-4, TrES-3, TrES-2 and WASP-3. For each 
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Fig. 20.— The transit times of WASP-3. P08: IPollacco et all (j20f)8T ): G08: lOibson et alJ (|200Sj ); 
EPOXI: this paper. Lower panel: An expanded view of the EPOXI transit times, with la errors 
of 23-51 seconds. 



system we provided an updated set of system parameters and orbital period, and placed upper 
limits on the secondary eclipse depth. For TrES-3, we see evidence of stellar variability over long 
timescales. We combined the EPOCh secondary eclipse upper limit for TrES-3 with previously 
published measurements and confirm that the data are best fit using an atmosphere model with 
no temperature inversion. For TrES-2, the EPOCh data weaken the cla imed trends of decreasing 
inclination and transit duration (jMislis &: Schmittll2009l : iMislis et al.ll2010l ). We have also performed 
a search for additional transiting planets in the EPOCh photometry for these systems, which we 
will present in a forthcoming paper. 
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Fig. 21. — The six EPOCh transits of WASP-3. The best fit model for the combined set of transits 
is plotted in each case. The scatter around the model is not typically larger in transit than out of 
transit, indicating that the in-transit residuals cannot be attributed to star spots. 
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